A MAP ON THE SPACE OF RATIONAL FUNCTIONS 



G. BOROS, J. LITTLE, V. MOLL, E. MOSTEIG, AND R. STANLEY 

Abstract. We describe dynamical properties of a map J defined on 
the space of rational functions. The fixed points of 5 are classified and 
the long time behavior of a subclass is described in terms of Eulerian 
polynomials. 



1. Introduction 

Let 9^ denote the space of rational functions with complex coefficients. 
The Taylor expansion at x = of i? G $K is written as 

(LI) Rix) = 

where n ^ —oo denotes the fact that the coefficients vanish for large negative 
n. We consider the map 5^ : ^ defined by 

(L2) d{R{x)) = Yl «2n+ix". 

The map ^ can be given explicitly by 



(L3) d{R{x)) 



2^ 



and it appeared in this form in the description of a general procedure for 
the exact integration of rational functions The splitting of an arbitrary 
function R into its even and odd parts R{x) = Re{x) + Ro{x) yields 

/•oo poo fOO 

(1.4) / R{x)dx = / Re{x)dx+ / Ro{x)dx. 

Jo Jo Jo 

The integral of the even part can be analyzed with the methods described 
in and the integral of the odd part can be transformed by x i-^ ^/x to 
produce 



poo poo poo 

(1.5) / R{x)dx = / Reix)dx + ^ / ^(Rix)) 

Jo Jo Jo 



dx. 
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Here we consider dynamical properties of the map : 9^ — > 9^. Section |21 
characterizes rational functions R for which the orbit 



(1.6) 



Orb(i?) := ^d^''\R{x)) : G n} 



ends at the fixed point G 9^. Section IHl describes the dynamics of a special 
class of functions with all their poles restricted to the unit circle. We estab- 
lish explicit formulas for the asymptotic behavior of their orbits, expressed 
in terms of Eulerian polynomials Am{x) defined by the generating function 



(1.7) 



1 — X 



y^Amix)—. 

z — / rn.' 
m=0 



1 — X exp[A(l — x)] ^ — 1 ' m 
The proof only employs the classical result (51, page 243: 



(1.^ 



Am{x) = {i-xr+'j2^ 



k=0 



Section m contains a description of all the fixed points of ^. 
The map ^ can be supplemented by 



(1.9) 



(B{R{x)) :-- 



for which similar results can be established. See |^ for details. These 
transformations can be written as 

= «'-'^^) and e(H(.)) = «(^-'^') +^^(^^^'^^) 

where L02 = —1- The extensions to higher degree will be considered in future 
work. For instance, in the case of degree 3 we consider the maps: 

1 



1i(ii(x)) := 



1 



3x1/3 
1 



3x2/3 



R{x^/^) + R{lo3X^/^) + R{ujIx^^'^) 

R[x^'^) + ulR{u^x^'^) + uj^Riujlx^'"^) 
R[x^'^) + u^Rioj^x^'"^) + ujIR{ojIx^'^) 



where = e^^^'/^. These maps correspond to subseries of the expansion of 
R taken along a fixed class modulo 3. 

Notation: For n G N, mi, • • • , m„ are odd integers. 

L = {mim2 • • • run)^^ and m* = mi + m2 + • • • + m„ — 2. 

For j £ and n G N 

71 

(1.10) Qn{x) = Hix^'^-l) 

k=l 
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2. The kernel of the iterates 

Even rational functions can be characterized as the elements of the kernel 
of ^. In this section we characterize those functions that vanish precisely 
after n applications of ^. 

The sets 

Kn := Ker^(") = ^Se^: ;?(")(S) = o} 

form a nested sequence of vector spaces. We now describe the class of 

functions J„ := Kn — Kn-i, i-c. those functions that vanish after precisely 
n applications of the map ^. In particular we show that J„ is not empty. 

The decomposition oi S e^R into its even and odd parts can be expressed 

as 

(2.1) S{x) = Si,iix^) + xS2,i{x^) 

where Si^i, 82,1 G $H. This decomposition applied to Si^i, 82,1 yields 

8ix) = 8i,2{x^) + X82,2{X^) + X''83,2{X^) + X=^S4,2(X^). 

Iterating this procedure produces a general decomposition. 

Lemma 2.1. Given n G N and 5 G 9^ there is a unique set of rational 
functions {-Sj,„ : < j < 2" — 1} such that 

(2.2) 8{x) = ^x^-i5,,„(x2"). 

Proof. Split the sum in 

(2.3) 8{x) = Yl 

fc3>— oo 

according to the residue of k modulo 2". □ 

We now show that the functions 8 e Jn are precisely those for which 
S2'\n = and S'2n-i „ 7^ 0. This generalizes the case n = 1 that states that 
d{S) = precisely when 8 is even. 



and 
(1.11) 
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Theorem 2.2. The rational functions that vanish after precisely n appli- 
cations of 5" are those of the form 

2"-l 

(2.4) S{x) = Y,x^-^S,,n{x^"), 

where Si^n, ■ ■ ■ , 'S'2"-2,n are arbitrary rational functions and 5'2"-i,ra 7^ 0. 
Proof. A direct calculation from ()2.2() shows that for 1 < A;, 

(2.5) d^'\S)ix) = 5^ x^-^S^.,- Jx2"-'=). 

i=i 

The statement now follows from 

(2.6) 5("-l)(5)(x) = 52n-l,Jx2)+x52n,„(x2) 

and 

(2.7) ^("^(5)(x) = 52n,„(x). 

□ 

Note. We now state the result of Theorem 12. 21 in the language of dynamical 
systems. Let f : X ^ X he a map on a set X and xq £ X he a fixed point 
of /, i.e. f{xo) = xq. We say that a sequence {xi, X2, • • ■ , x„} of elements 
of X is a prefixed sequence of length n attached to xq if Xj+i = f{xi) and 
f{xn) = Xq. Theorem 12.21 states that ^ admits prefixed sequences attached 
to of arbitrary length. 

3. The dynamics of a special class 

The asymptotic behavior of 5^ can be described in complete detail for 
rational functions in the class 

f Fix) 

(3.1) (mi, . . . , nin) = { ^ , X : -P is a Laurent polynomial 

where mi, • • • , m^ are odd positive integers and Qn{x) is defined in H1.1U() . A 

Laurent polynomial is a rational function of the form a--kX~^ +a-k+ix^^^'^ + 
• • • + Oj-ix^^"^ + OjX^ with k, j G N. 

The case n = 1 has been described in f?. The results are expressed in 
terms of the function 



(3.2) 7„(j) = m 



i(m-l)(j-l) 



(m — 1 + j)/2 if j is even 
{j-l)/2 ifjisodd. 
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Theorem 3.1. Let m be an odd positive integer. Then 

1) For j € N 

Thus, the study of the dynamics of ^ on (^i is reduced to that of ■jm on Z. 

2) The iterates ^'ym\j) : p = 0, 1, 2, • • • | reach the set 

(3.4) 21^ := {0, 1, 2, ••• ,m-2} 

or the fixed points — 1 and m — 1 in a finite number of steps. Moreover, 
is invariant under the action of jm- This action partitions 21^ into orbits. 

3) The inverse of the restriction of 7^ to 21^ is given by 

\ 2k + 1 — m it [m — I) / 2 < k < m — 2. 



Note 3.2. The exphcit form of 6m permits the exphcit computation of the 
orbits of 7m on the invariant set 2lm- For example, if m is prime then every 
orbit of 7m is of length Ord(2;m). In particular, there is a single orbit if 
and only if 2 is a primitive root modulo m. 



The iterates 7m (j) can be characterized by the congruence (|3.6|) below. 
This will be used in the determination of the limiting behavior of the iterates 
of ^ below. The proof of this congruence and the numerical and symbolic 
evidence of the asymptotic behavior of the iterates of ^ on Rj^n{x) were part 
of a SIMU 2002 project. Details will appear in (4,. 



Lemma 3.3. Let m be an odd positive number, < j < m, and p G N. 
The unique solution of 

(3.6) 2P{x + 1) = j + 1 mod m 

in < X < m is given by x = 7m^(j)- 

Proof. The proof is by induction on p. Note first that the solution of the 
congruence (|3.6j) is unique because gcd(m, 2^) = 1. The base case {p = 0) is 

X + 1 = J + 1 mod m 
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with unique solution j = 7m^(j)- To complete the inductive step observe 
that 



m 



m 



Iraij) 



-2^m (^^P)(j)-l)+2^'(7^^)(i) + l) 



= j + 1 mod m. 



□ 



The map has a very rich dynamical structure, even in the case n = 1. 

Theorem 3.4. Let r G N. Then ^ has at least one periodic orbit of length 
r. 

Proof. For r G N define m = T' — 1. Then the orbit of l/(a;"* — 1) under 5^ 
is 



X 



xm _l x'^ -1 



□ 



We now consider the properties of the class for n > 1. 



Lemma 3.5. The class £„ is invariant under the action of S^. 

Proof. We show that ^{Rn,j{x)) G £„• The linearity of ^ yields the result. 

Introduce the notation 

(3.7) 5^ ^a,(mii+mi2+-+mi,)/2 

where the sum is taken over all subsets of {mi, • • • ,m„} containing I ele- 
ments, the empty sum giving S'o = 1. 
Now 



(3.8) 

and 
(3.9) 



n n 



k=l 



n n 

= (-l)n + ^(-l)n-^5,. 
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SO that 

=. X S 

The sum S contains all the terms Si with / of the opposite parity of j so 
that x'^-'^"'^-*/^ X S" is a polynomial. □ 

We now consider the orbit of a general rational function R in the class 
The case n = 2 illustrates the general situation. A direct calculation 
shows that 



/ \ ^(i-l+mi)/2 _^ 2,(j-l+r7i2)/2 

3^ I -, r-, r I = ^ if 7 is even 

\{x'^^ - l){x'^^ - I) J - l)(x™2 - 1) 



and 



a;0-l)/2 _^ rj.(j-l+mi+m2)/2 

if j is odd. 



- l)(a;™2 - 1) y - l){x'^2 _ 1) 

Thus 5^ preserves the denominator of R (as was shown in Lemma l3.5() and 
each monomial of R yields two monomials. We now show that the exponents 
of these monomials are always bounded. 

n 

Proposition 3.6. Let m* = ^ — 2 and define 

k=l 

r p 

21 := < — - G : P is a polynomial with deg(P) < m* 

Then 21 is invariant under 5^, and every orbit starting at R ^ <tn reaches it 
in a finite number of steps. 

Proof. Let R = P/Qn S 'tn and write 5'(i?) = Pi/Qn- The largest exponent 
in Pi, call it j, appears from the sum Sn- The inequality j > m* implies 
^(j — 1 + mi + • • • + m„) < j. The case j < is similar. □ 

Now we consider the asymptotic behavior of the iterates of ^ starting 
at i? € This is expressed in terms of the Eulerian polynomials Am{x) 
defined in (|1.7|) . The discussion is divided into two cases depending whether 
the number 

(3.10) d := gcd(mi, m2, • • • ,"?-n) 
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is 1 or not. In Theorem 13 . 71 we prove that if d = 1 then 

/ N 2("~^)P A 1 (r) 

(3.11) d^^HRix)) ~ ^- ^p^^ 

^ ^ u V V (n - l)!mi • • -m^ (1 - x)" 

The case d > 1 is described in Theorem 13.81 We prove that the sequence of 
iterates of ^ apphed to the function x^ /Qn{x) taken along a fixed residue 
class (defined in (|3.1<S|) ) has an asymptotic behavior as in the case d = 1, 
with limit points expressed in terms of Eulerian polynomials. 

This suggests the existence of an arbitrary number of limit points, but we 
have not ruled out the possibility that all these could coincide. 



The proof employs the observation that if 
(3.12) R{x) = Y,f{k) 



x"^ 



is the expansion of i?, then 

(3.13) ^^P\R){x) = /(2^(fc + 1) - 1)^'- 

k 

Theorem 3.7. Let i? € 21 and suppose d= 1. Then 

(3.14) ita?!^ = , ^f;'"' 

^ ' P^oo 2("-i)P (1 - x)*^ (n - 1)! 

where L = l/{mi ■ ■ ■ m„). 

Proof. Since d = 1, the rational function 

P{x) 



(3.15) R{x) = , 

has a pole of order n at x = 1 and poles of order less than n at the other 
zeros of Qn{x\ all of which are roots of unity. Thus the partial fraction 
expansion of R has the form 

(3.16) R{x) = ^^^-^ + G(x) 

where the term G{x) contains all the terms of order lower than n. Hence 

(3.17) R{x) = lJ2C^^'^~^)x'' + G{x) 

k>0 ^ ^ 

where the coefficient of x^ in G{x) is 0{k"'~'^). Thus 



^(P\R){X) = Lj;(^ + ''^^_+^2'-')x'=+^?(^)(G)(x) 



[n — 1)1 ^-^ 
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where we have used 

'n + 2Pk + 2P-2\ 2("-i)P 



n — 1 J (n — 1)! 

as p ^ oo. The result now follows from (|1.8|) . □ 

The analysis of the case d > 1 involves the function 

with d as in (|3.1U|) . The function p appears in the dynamics of 'jm- for 
< j < m — 2 the length of the orbit containing j is p{j)- See [2] for details. 

Introduce the notation 

(3.19) Ei := {j eZ: there exists r G N : tJ'Hj) = i } 
for the backward orbit of and let 

(3.20) A := E^iUEd-i 

be the integers that eventually are mapped to the fixed points of 7^. 

Theorem 3.8. Let n e N. 

a) If j G ^ then 

lim 



p^oo 2("-i)P (n - 1)! (x'^ - !)"■ 

b) If j ^ A choose r* such that 

r :=7f)(j)G{0, 1, 
Then for any q £ {0, I, ■ ■ ■ , p{f) - 1} 

j.^ g(PP(^"')+g)(fi,,„(x)) 



where 7 = 'y'f\j*)- 

Proof. Since gcd(mi, m2, • " " I'^^n) = the function Rj,n{x) has a pole of 
order n at the d-th roots of unity and a pole of strictly lower order at all 
other zeros of Qn{x). Thus its partial fraction decomposition has the form 

R-i n(x) = -, h 7 ^^^-T !-••• + -, ^ "^'"^ h lower order terms 

^'^ ^ (x-l)" (x-wi)" (x-Wd-i)" 
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where ui = exp(27ri//(i), 1 < I < d — 1. The coefficients Cij are given by 



Thus 

(3.21) Rax) = i:7-^^ + G(x), 

1=0 ^^'^ 

where the error term G{x) includes the polar parts of all poles of order less 
than n. Then 



oo d-1 / . J i\ 
k=0 1=0 ^ ^ 



Therefore 



k=0 \i=o 

Now 

d-l 



if d does not divide q 
d ii d divides q, 



1=0 

so the only values of k that contribute to the sum are those for which 

j + 1 = 2'^(A; + 1) modd. 

The discussion is divided into two cases depending on whether some iterate 
of j reaches one of the fixed points or not. 

Case 1: Suppose j € A and let r* G N be such that 7^^ \j) = —1. Then 
Lemma l3 . 31 shows that k = ^^J^^ \j) + Nd = Nd — 1. Thus 

5(^+'^*)(i?,,„(x)) = 

N=l ^ n 1 J 

In the limit as r — > 00 the binomial coefficient is asymptotic to 

2{r+r*){n-l) j\j-n-l^n-l 

(n-1)! ■ 
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The number r* is fixed, so we have 

^' ' r^oo 2(^+'-*)('^-i) (n - 1)! (1 - x'^)'^ 

as stated. 

A similar argument shows that the same result is true if j lies on the 
backward orbit of the second fixed point. 

Case 2: Now assume j ^ A. Lemma 13.31 shows that 
Thus 

N=o ^ n - 1 J 

In the limit as r — > 00, the binomial coefficient is asymptotic to 

(n-1)! 

where 71 = 7'-^"'"^*''(j)- Thus 

^ ■ N=0 

as r — > cxD. Now 

00 n—l 



"^f-" Af=0 1=0 ^ ^ 



Af=0 Af=0 1=0 

n-1 



1=0 

As a function of r, 

7i=7^^+''*^(j)=7^^Hj*) 
has period p{j*)- Write r = pp{j*) + g with < g < — 1 and replace 

7^^(r)by7?^(r)to obtain 

2(pp0*)+9+'-*)(™-i) ~ 
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To conclude the proof, observe that as q runs over the set of residues modulo 
pij*), so does q + r*. □ 

4. The fixed points of 5" 

A formal power series argument shows that any rational function R fixed 
by ^ must have an expansion of the form 



n=0 



oo 



(4.1) xR{x) = c + ^a„(/p(x2"+i), 
where 

(4.2) ip{x) = ^x^* =x + x^ + x'^ + x*H . 

n=0 

In particular, R has at most a simple pole at the origin and if 

n>-l 

is such a fixed point, then f{2n) = f{n) for n > 0. Thus the problem of 
finding fixed points of ^ is reduced to finding sequences {a„} for which 1)4. 1() 
is a rational function. 

The class of functions discussed in Section |31 yields examples of fixed 
points. Let m be an odd positive integer. Then m — 1 is fixed by 7m, so 
Ri,m-i{x) = x"^~^/{x"^ — 1) is fixed by ^. This example can be obtained 
by a different approach. First observe that if ^{R) = R and r is any odd 
positive integer, then the function 

(4.3) ^r{Rix)) = x'-^R{x'') 

is also fixed by ^. The function g{x) = l/(x — 1) is fixed by ^, so that 
Ri,m-i{x) = ^rn{g{x)) is also fixed. 

The description of all the fixed points of ^ requires the notion of cyclo- 
tomic cosets: given n, r G N with r odd and < n < r — 1, the set 

(4.4) Cr,n = mod r : s G Z} 
is the 2-cyclotomic coset of n mod v. Observe that Cr,n 

is a finite set. With 

A a fixed primitive r-th root of unity, define 

(4-5) frA-) = E r^A^^- 

The partial fraction decomposition of the fixed point ^rnisix)) can be 
decomposed into a sum of rational functions each fixed by 5- For example. 
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consider ^Y{g{x)) = jix^ — 1) and let A = exp(27ri/7) be a primitive 7-th 
root of unity. Then 

6 



(4.6) 7*B7(5(x)) = Y.\ 



X^x 

k=0 



A A^ 

+ I Z ^ + Z + Z I + 



1-x \1-Xx l-X'^x l-A^a 

A3 A5 A6 
+ + 



- X^x 1-X^x l-A^a;^ 

where each of the sets of terms grouped together is a rational function fixed 
by ^. In the notation introduced above, this decomposition is 

(4.7) 7^7{g{x)) = /7,o(a;) + /7,i(x) + /7,3(x). 

We now classify all the fixed points of 5- 



Theorem 4.1. A rational function is fixed by 5 if and only if it is a linear 
combination of ^ and the functions fr,n{x) for r odd and < n < r — 1. 

Proof. The identity 

A 



l-Xx J 1-X'^x 

shows that ^ fixes the fr,n because the squaring map S : z i-^ z"^ permutes 
the values A™ for m G Cr,n- 

We first establish the converse under the assumption that the poles of R 
are simple. The final step of the proof consists of checking that this condi- 
tion holds for any fixed point of 3^. 

Let be a rational function, with simple poles, that is fixed by The 
partial fraction decomposition of R is 

which is unique up to order. Apply ^ to produce 

c ^ ajX] 



]=l 3 



R{x) = -+y- 

X ^ 1 

The uniqueness of (|4.8|) shows that 

(4.9) {Ai,---,Aj} = {Xl---,X'j}, 

that is, the set A = {Ai,-- - ,Aj} is permuted by the squaring map. We 
conclude that the set {A| : n G N} is a finite set (a subset of A) and so 
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every is an r^-th root of unity, for some odd positive integer r^. 

Now group terms in the sum ()4.8() according to the orbits of the squar- 
ing map S on the set A. The coefficient aj in (|4.8j) must be constant 
along each orbit, and moreover, the orbit of under S is precisely the 
set {A™ : m G Cn,rk} for some n. Therefore R{x) can be decomposed as a 
linear combination of the required form. 

The next result concludes the proof of the theorem. 

Proposition 4.2. Let R{x) be a rational function that can be expressed in 
the form 

oo 

(4.10) xR{x) = c + ^a„99(x2"+i) 

n=0 

where (p is given in (|4.2j) . Then the poles of R must be simple. 

The following result, used in the proof of Proposition lT^ is demonstrated 
in [HI, page 202. 

Lemma 4.3. Let qi,q2, ■ ■ ■ ,Qd be a fixed sequence of complex numbers, 
d > 1, and ^ 0. The following conditions on a function / : N — > C are 
equivalent: 

(1) Eti>o/(^)^" = ^ where, Q{x) = 1 + qix + q2X^ + q^x^ ^ hqd.x'^- 

(2) For > 0, 

k 

/(n) = ^P.(n)Ar, 

i=l 

where l + qix + q2X^ + qsx^ H h g^x'' = ni=i(l " the Ai's 

are distinct, and Pi{n) is a polynomial in n of degree less than dj. 

Proof of Prop. 14. 2t Assume P{x),Q{x) are relatively prime, and that f{n) 
is the generating function for P{x)/Q{x) written in the form promised by 
the lemma above (for n ^ 0). Since f{n) = /(2n), 

k k 

Q(x) = n(i - A.x)* = na - A?xr% 

i=l 1=1 

so 

{Ai, . . . , Afc} = {Xf, ... , Afc}. 

As in the proof of Theorem l4.11 we conclude that each Aj is a primitive rj-th 
root of unity for some positive integers ri, . . . , r^. 
Let M = lcm(ri, . . . ^rj.) and for a G N, define 

Ra = {fn G N : m = a mod M} 
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and fa 
Then 



to be the restriction of the function / : N ^ C to the set Ra- 

Ra 



faia + jM) = 5^ P.(a + jM)Xl+^'^ = ^ + jM)Af , 



so each fa has a representation as a polynomial in the variable j since A" 
is constant on the set Ra- We denote the natural extension of this map to 
an element of C[j] by Fa- Note that the restriction of Fa to N will not be 
/ in general. Our goal is to prove that each Fa is a constant function, with 
corresponding constant denoted by Ca- Once this is shown, we have 

Q{x) Z^'^^Z^^ Z^i_^M' 

^ ^ ' a=l j=0 a=l 

SO P{x)/Q{x) is a rational function with only simple poles, as desired. 

It remains to show that each polynomial map : C ^ C is a constant 
function. For each positive integer n, define Sn = {2*n : t G N}. We say 
that a has an infinite cross-section if Ra H Sn is an infinite set for some 
n G N. We proceed by considering two cases, depending on whether a has 
an infinite cross-section or not. 



Case 1: Suppose a has an infinite cross-section, i.e., Ra H Sn is an infinite 
set. Since f{s) = /(2s) for all s G N , is constant on Ra n 5„. Since 
Ra n Sn is an infinite set, Fa is a constant polynomial. 



Case 2: Suppose a does not have an infinite cross-section, i.e., Ra H Sn 
is finite for all positive integers n. Then Ra H Sn must be nonempty for 
infinitely many values of n. Since there are only finitely many distinct sets 
of the form i?;,) it follows that for each S^ there exists 6 G N such that 
-Rfo n Sn is infinite. Moreover, since there are only finitely many choices for 
-Rfc, there is at least one 6 G N such that there exist infinitely many values of 
n where Ra H Sn is nonempty and R^ H Sn is infinite. Since b has an infinite 
cross-section, an application of Case 1 demonstrates that the restriction of 
/ to Ri) is the constant function c^. Since / is constant on each Sn, the 
restriction of / to Sn is the constant C},. Thus Fa achieves the value cj, 
infinitely many times, and so Fa must be a constant polynomial. □ 

The proof of Theorem 14. II is complete. 

Note 4.4. For each fixed point Rf we construct the rational function 

2"-l 

(4.11) R{x) = ^x^-'Rj^n{x^")-Rf{x), 
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where Ri^n,''' ,R2"-2,n are arbitrary rational functions and ii2"-i,n 7^ 0. 
Theorem 12.21 shows that R is the general form of a prefixed sequence of 
length n attached to Rf. 
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